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Abstract Fe-based superconductors (FeSC) exhibit all the properties of systems 
that allow the formation of a superconducting phase with oscillating order param- 
eter, called the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. By the analysis 
^ , of the Cooper pair susceptibility in two-band FeSC, such systems are shown to 

O ■ support the existence of a FFLO phase, regardless of the exhibited order param- 

eter symmetry We also show the state with nonzero Cooper pair momentum, in 



?-H ■ superconducting FeSC with ^ cos{kx) • cos{ky) symmetry, to be the ground state 



of the system in a certain parameter range. 
Keywords FFLO • pnictides 
PACS 74.20.-z • 74.70.Xa • 74.81.-g 



C ■ 1 Introduction 

o 
o 



At low temperatures the orbital pair breaking effects are smaller in magnitude 
than the Pauli paramagnetic effect, so that superconductivity survives up to the 
Pauli limit - a phase with oscillating order parameter (called the Fulde-Ferrell- 
Larkin-Ovchinnikov phase or FFLO in short) [T] can be more stable than a phase 
with a constant order parameter (the Bardeen-Cooper-Schrieffer phase, or BCS in 

^y^ ■ short). In this case. Cooper pairs may be formed with non-zero total momentum 

^S| . between Zeeman-split parts of the Fermi surface. 

Properties of this phase have been usually evaluated in tight-binding models of 

^^ ', one-band systems. [2j However the latest experimental [3j and theoretical [4J works 

^T) ■ suggest we can expect the existence of the FFLO phase in multi-band Fe-based 

superconductors (FeSC). It follows from the fact that they possess properties close 
to heavy fermions systems, [5j for which strong experimental evidence suggest the 
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existence of said phase. [6 Both kinds of systems are mult i- layered, clean and have 
a relatively high Maki parameter. 

In this paper, making use of the Cooper pair susceptibility and the the mini- 
mization of free energy of the system, we discuss the possible appearance of the 
FFLO phase in pnictides. In Section [2] we describe the selected model of FeSC, 
in Section [3] we present our methods. In Section |4] we illustrate and discuss our 
numerical results. We summarize the results in Section O 



2 Theoretical model 

The FeSC system is described using a two-orbital per site model, with hybridiza- 
tion between the dxz and dyz orbitals. We adopt the band structure proposed in 
Ref. [7] and assume that the external magnetic field is parallel to the plane. The 
Hamiltonian of the system in momentum space takes the following form: 

^0 = ^ X^(^fc ^ - (m + ^h)6ai3)clj^^c^ka (1) 

ka af3 

Tf^ = —2 (ti cos{kx) + t2 cos{ky)) — Ats cos{kx) cos{ky), 
Tp, = —2 (t2 cos{kx) + ti cos{ky)) — Ats cos{kx) cos{ky), 
Tk^ = Tk^ = -4t4sin(/ca,)sin(^), 

where c^^j^^ (c^^ka) is the creation (annihilation) operator of a particle with mo- 
mentum k and spin a in the orbital a. T^^ is the kinetic energy term of a 
particle with momentum k changing the orbital from /3 to a, /x is the chemical 
potential and h is the external magnetic field. The hoppings have magnitudes: 
(ti, t2, t3, t4) = (—1.0, 1.3, —0.85, —0.85), in units of |ti|. At half-filling, a configura- 
tions with two electrons per site requires /i = 1.54|ti|. Our choice of the parameter 
set is motivated by the fact that it reproduces the same Fermi surface structure 
as the local-density approximation calculations of band structure. ^Sj 
By diagonalizing the above Hamiltonian, one obtains 



^0 — 2_^E^kad^j^^d^ka (2) 

ekcr 

with eigenvalues E^ka = ^ek ~ (m + cr/i), where: 




E±,u = ^-^^^ ± \l { ^-^^X+ [n-^f, (3) 



c/^^^ is a new fermion quasi-particle operator in the band e = =b. In this case we 
have two Fermi surfaces (Fig. HJa) - giving an electron-like band [e = +) and 
hole-like band [e = —). 
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Fig. 1 (Color online) The vector <5 defines the pairing between sites i and i -\- 5 for different 
symmetries of the order parameter. Colors and symbols correspond to the sign of the order 
parameter for a given direction in real space. For s-wave symmetry the pairing is between two 
electrons on the same site of the lattice, while for other symmetries it is between two other 
sites (nearest neighbors or next nearest neighbors). In contrast to d type symmetries, s type 
symmetries do not change sign depending on the direction. 



3 Methods 



We introduce a superconducting pairing between the long-lived quasi-particles in 
bands e = ±. [9] To determine the possibility of formation of the FFLO phase, we 
turn our attention to the static Cooper pairs susceptibility: 






£J' 



(4) 



where ((. . .)Y is the retarded Green's function and A^i = ^ . 'd{j — i)d^i^d^j^ is 
the OP in band s. The operator d^i^j. in real space corresponds to the operator 
c/gfccr in momentum space. The Factor '^{j — i) defines the OP symmetries (Fig.[T]) 
- for example for d^2_y2-wave pairing, 'd{5) is equal to 1 ( — 1) for ^ = ±x (±^) and 
zero otherwise. In momentum space: 



(g) = Wm. --^r]{-k - q)T]{l)Gs{k,l,q,u), 



(5) 



hi 



Gsik,l,q,uj) = {{dsktde-k-qi\dl_i_q^dli^)Y 



^kl- 



f{-Esk^) - f{E^_j^_q^) 



^£fct 



-h-qi 



(6) 



where ri{k) is the structure factor: 



7?(fc) 



1 for s-wave 

2 (cos(A;a;) + cos(ky)) for S2^2^^2-wave, 

4 cos{kx) cos(ky) for Sj,2y2 (s±)-wave, 
2 (cos(A:a;) — cos{ky)) for da,2_^2-wave, 

4 sin( /ca; ) sm(ky ) for 



(7) 
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corresponding to the type of symmetry of the OP. 

We investigate the tendency to form the FFLO phase in the system, using 
the static Cooper pairs susceptibihty Xe (q)- I^i magnetic fields of the order of the 
Pauh Hmit, when the critical FFLO field {he ) is bigger than the corresponding 

BCS field (hf^^), the FFLO phase is favored. In such case, the divergence of 
this function for some q y^ may imply a second-order transition to the FFLO 
state of corresponding symmetry from the normal phase. [TO^ The location of the 
maximum of the response function Xe (q) matches the preferred momentum of the 
Cooper pairs in the system described by the Hamiltonian ([2]) in magnetic field h. 
This method allows to establish the propensity to form the superconducting phase 
(with non-zero momentum of the Cooper pairs) without specifying the mechanisms 
responsible for the ordered phases with given symmetry. Additionally we obtain 
the change in the pair susceptibilities Sx^iq) = Xs^iQ) — xf^o) due to the external 
magnetic field {xeio) with the field, Xeio) without respectively). 

It should be noted that the divergence of the Cooper-pair susceptibility is 
neither a sufficient condition nor evidence for the transition to the FFLO state. In 
order for this to happen the system energy f2{q) should attain its minimum at a 
nonzero Cooper pair momentum q in a magnetic field h > he , equivalent to the 
condition hc^^^ > hf^^. To check this, we effectively describe superconductivity 
in the FFLO phase by the Hamiltonian: 

HSC = ^ {^ekdsktde-k^qel + ^-C-) , (8) 

ek 

where A^j^ = Asr]{k) is the amplitude of the OP for Cooper pairs with total 
momentum qs (in band e with symmetry described by r]{k)). As we see, in the 
operator basis d^j^^j. the total Hamiltonian H = Hq + Hgc formally describes a 
system with two independent bands. Using the Bogoliubov transformation we can 
ffnd the eigenvalues of H: 



Xf, = ^-^-^-■-^+'^^ ±./( ^-H + ^---''+-^^ )%|z^,,|2. (9) 



The free energy is given by: 



^ = -^^ E E 1^ (1 + exp(-/3A,\)) + ^ (Sefc; - ^^) , (10) 



aG± ek ek 



where Ve is the interaction intensity in band s. The global ground state for ffxed 
h and T is found by minimizing the free energy w.r.t. the OPs and q. 



4 Numerical results and discussion 

Numerical calculations were carried out for a square lattice Nx x Ny = 600 x 600 
with periodic boundary conditions, for kT = 10~^|ti|. As a ffrst step the static 



Cooper pairs susceptibility Xsiq) was calculated in magnetic ffeld h = 0.25 |ti|. 
Then the free energy i7(qr) of the superconducting system was evaluated for mag- 
netic ffelds near the Pauli limit hp c^ 0.25 |ti|. 
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Fig. 2 (Color online) The static Cooper pairs susceptibility xf'io) ii^ the presence of the 
external magnetic field h = 0.25|ti| and kT = 10~^|ti| for different symmetries. 
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Fig. 3 (Color online) Change in the static Cooper pairs susceptibility Sx^{q) (for data pre- 
sented in Fig. [2} 
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The static Cooper pairs susceptibility. Assuming different symmetries 77 (fc) of the 
superconducting OP in bands £ = it, we characterized the Cooper pair suscepti- 
bility - Fig. [2j For every OP symmetry, in the band e = — the static Cooper pairs 
susceptibility X- (9) takes its maximum for g 7^ 0. Conversely in the band e = +, 
with Sx2_^y2 and dj,2y2 symmetry of the OP, there is a strong tendency to form a 
BCS phase (maximum x+ (q) for q = 0). When he > he this can be a sign 

of the existence of the FFLO phase in the band s = —, while the band e = + is in 
the normal state. Numerical data for both d-wave type symmetries in band s = — 
is less clear cut, as the maximum x^q) is only slightly greater than x(0)- 

There is a clear preference in case of £ = — for much smaller momenta than in 
band e = +, due to the relative width of the bands. Cooper pair momenta depend 
on the split in the Fermi surface, caused by the external magnetic field, which is 
larger for the broader band e = +. Additionally the presence of a magnetic field 
causes a dampening in each case of the response function near zero momentum 
(Fig. [3]- in blue). Nonetheless larger momenta are unaffected and increasing (in 
red). 

The behaviour of the response function Xe (q) shows that multi-band systems 
have the characteristics typical of one-band systems - Cooper pairs in the FFLO 
phase possessing momentum along the principal directions of the system are pre- 
ferred, [11] - for example in directions [±1,0] and [±1,1] for s-wave and dj.2_y2-wave 
symmetry respectively in band £ = + (Fig. [2]). 
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Fig. 4 (Color online) Detailed study of the minimal two-band model describing iron-base 
superconductors with s^2y2{s±)-wave symmetry proposed by Ref. [7 . (Panel a) Fermi surfaces 
(solid line) for fi = 1.54|ti |. The background color describes the sign of the OP (red for ri{k) > 
and blue for r]{k) < 0). (Panel b) The free energy i7(g) in the two bands £ = ±, for different 
values of the Cooper pair momentum q, showing the location of the minima and indicating 
the existence of different phases. Results for h = 0.25|fi| and temperature kT = 10~^|fi|. 



Minimization of free energy. Theoretical results indicate the presence of 3^2 y2 ^ 
cos{kx) ■ cos{ky) (also called s±) pairing symmetry in FeSC . [T^ In this case the 
OP exhibits a sign reversal between the hole pockets and electron pockets (Fig. 
Sla). Taking this into account, in this paragraph only consider 53,2^2 symmetry. Vs 
was taken such that the Pauli limit was of the order hp :^ 0.25 |ti| (V+ = — 0.74|ti| 
and V- = -1.56|ti|). 

The study of the free energy f2{q) for the BCS state {q = 0) w.r.t. magnetic 
fields h c:^ hp, showed that phase transitions in both bands are first-order for 

2 and dj,2y2 which are second-order. Only the 



S^2 



x^-\-y 



all symmetries, except for 
minimization of f^{q) w.r.t. qr, allows to check whether the system exhibits a FFLO 
phase. Varying q G FBZ in case of 53,2^2 pairing showed that the band £ = -\- 
undergoes a transition from BCS to the normal state and the band e = — from 
BCS to FFLO state (Fig. Hlb). Further increasing the magnetic field, the FFLO 
phase persists in e = — . It should be pointed out that in this band exist four 
equivalent Cooper pair momenta (±q, 0) and (0,±^), in agreement with the static 
Cooper pairs susceptibility results, and also with previous works. [TT] Moreover, 
it is reasonable to expect that the phase with an OP given by the superposition of 
plane waves with said momenta would be energetically favored by the system. [T3] 



5 Summary 



FeSC exhibit many characteristic features of systems in which we can expect the 
existence of the FFLO phase. Using a minimal two-band model for FeSC, we 
conducted a numerical study of FFLO phase in multi-band systems. The static 
Cooper pair susceptibility suggests that we can expect the system to prefer the 
state with nonzero Cooper pair momenta (the FFLO phase) regardless of the OP 
symmetry, when hc^^^ > hf^^. Moreover, the ground state of the system with 
Sx2y2 ^ cos{kx) • cos{ky) symmetry OP, can be the state with nonzero Cooper pair 
momentum for magnetic fields near the Pauli limit. 
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